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1. Introduction. 

To this author, limit theorems in probability theory have always been 
irresistibly attractive, as well as centrally important. As Kolmogorov put it 
in his introduction to Bernoulli's book, "The cognitive value of probability 
theory lies in the establishment of strict regularities resulting from the com- 
bined effects of mass random phenomena. The very notion of mathematical 
probability would have been fruitless if it were not realized as the frequency 
of a certain result under repeated experimentation. That is why the works 
of Pascal and Fermat can be viewed as only the prehistory of probability, 
while its true history begins with J. Bernoulli's law of large numbers" ([22]; 
[37], 874-5). The weak law of large numbers gives us a limiting frequency in 
the presence of a mean; its refinement, the central limit theorem, gives us a 
limiting distribution in the further presence of a variance. Both phenomena 
can occur more generally. The necessary and sufficient condition is known in 
each case, and in each case involves regular variation (for details see e.g. [4]). 
Thus probabilists are led inescapably towards regular variation (for which see 
e.g. [5], particularly Chapter 8). Furthermore, the role of the characteristic 
function, or Fourier transform, in translating averaging, or addition of inde- 
pendent random variables, into multiplication rather than convolution leads 
one to work with transforms rather than with distributions. But one needs 
eventually to translate the information thus obtained from transform to dis- 
tribution language. This is the role of Tauberian theory. For a magisterial 
treatment of the whole subject, see [24]; for regular variation, see Chapter 
IV there, or [5] Chapter 4. 

For background, we mention here what Tauberian theorems are. When 
handling limits, or asymptotic properties, we often need to pass from an ob- 
ject to some smoothed version of it. This is usually easy, and is an Abelian 
theorem (prototypes: Abel's continuity theorem for power series, or the impli- 
cation from s n —>■ c to a n — > c, with o n the arithmetic mean (si + . ■ . + s n )/n). 
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The converse is false in general, but may become true under some additional 
condition. Such a result is a Tauberian theorem, and the additional con- 
dition (the 'signature' of the Tauberian theorem) is a Tauberian condition. 
Sometimes one may have a hypothesis on both object and smoothed version, 
and a conclusion on the object, with no additional, or Tauberian, theorem. 
Such a result is a Mercerian theorem (prototype: Mercer's theorem, that 
(s n + °n)/2 -> c implies s n ^ c). 

Our starting point here is §4.12 of [5], on Tauberian theorems of exponen- 
tial type. There the treatment closely followed that of Kasahara [19] (leading 
to three Tauberian theorems, named after de Bruijn, Kohlbecker and Kasa- 
hara). At that time, neither Yuji Kasahara nor I realized the link between 
this area and large deviations. This link has been subsequently elaborated, 
by Kasahara and Kosugi; see [25], [26], [20], [21]. It involves convex analysis 
and the Gibbs variational principle; for background, see [40], [17], [11], [10]. 
The theme is exemplified by the relation displayed on the cover of [10]: 



(J is the rate function). 

For later use, we record here our notation for the Fenchel dual (or Legendre- 
Fenchel dual, or conjugate) of a convex function /: 



([36] Sections 7, 12). Then /* is convex, proper (< +oo somewhere, > 
— oo everywhere) iff / is proper; the closure (greatest lower semi-continuous 
minorant if / is nowhere — oo, identically — oo otherwise) cl / has the same 
dual as /, /* is closed (equal to its closure), and /** = cl /. 

We shall only be concerned here with the case when all moments fi n : = 
i?[X n ] of X are finite (in this case the moment-generating function (mgf) 
M(t) := E[e tx ] = Y^=o^nt n /n\ is defined, though its radius of convergence 
may be zero). For X non- negative (which we may suppose here without 
loss) the non-integer moments are then finite also, and one has the option 
of working with the Mellin transform (or Mellin-Stieltjes transform) T{s) : = 
E[X S ]. (Of course, we may pass between the mgf and the Mellin transform by 
an exponential, or logarithmic, change of random variable, but this may not 
be probabilistically natural.) Just as the cf is the natural tool for handling 
sums of independent random variables, the Mellin transform is the natural 




f*(x) := sup y {xy - f(y)} 
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tool for handling products. For a systematic treatment, see [41]; for some 
applications, see Bingham and Teugels [6]. We return to Mellin transforms 
in Section 4. 

2. Entire characteristic function 

We know from the work of Paley and Wiener [34] that the proper setting 
for the Fourier transform is the complex plane (or domain, in their termi- 
nology). Thus the traditional 'i' in the definition <j>(t) :— J e ltx dF(x) of the 
characteristic function (cf ) of a probability distribution F is merely conven- 
tional, and may be omitted if we intend to let the argument t take complex 
values. One then needs to distinguish cases according to the nature of the 
cf, or mgf, as a function of a complex variable. It may be entire, analytic but 
not entire or non-analytic, according as the radius of convergence is infinite, 
positive and finite, or zero. We turn now to the entire case, for which see 
e.g. [29], and [5] Section 7.6 (Levin-Pfluger theory: functions of completely 
regular growth). For the analytic case, see [29] Ch. 7, [30] Ch. 7. For the 
non-analytic case, see Section 3 below. 

An entire cf (of a positive random variable X, with distribution function 
F) has order p e [l,oo) ([29], Th. 7.1.3). If for large t \ogE{e tx )/t p has 
liminf a and limsup b with < a < b < oo, we say that X has very regular 
growth ([38] II. 7). Davies [9] shows that in this case one can find upper and 
lower bounds on 

-\ogP{X > x)/x p/{p ~ 1 \ 

in terms of a and b. If the limit exists, and a = b = a, say, then we say 
that X has completely regular growth. The Levin-Pfluger theory of [28], [5] 
Section 7.6 then applies, and — logP(X > x)/x p ^ p ~ 1 ^ is regularly varying 
with index p — p/(p — 1). 

This last result is extended to the case of general regular variation by 
Geluk [14] . Writing R p for the class of functions regularly varying with index 
p, he links 

lim inf — log P(X > x)/ <f>(x) = 1 

for (ft & R p with 

lim sup [E(X n )] l ' n ij(n)/n = q 1/q e- 1/p , 

x— >oo 

where p > 1, q > 1 are conjugate (in the usual sense of L p spaces) and tp is 
increasing and in R±/ q . Then 

ip~*(x) ~ <j)*(x), 
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where one has the inverse function of ip on the left and the Fenchel dual of 
on the right, so both sides are in R q (see e.g. [5] Section 1.8.4). 
Example: The Mittag-Leffler distribution. 

For a G (0,1), the Mittag-Leffler law with index (or parameter) a has 
mgf 

oo 

ML a (t) := 5>7r(l + na); 

n=0 

for uses in probability theory, see e.g. [3], [12] XIII. 8, [5] Section 8.11. Thus 

= n\/r(l + no) = r(n)/(aT(na)), 
and by Stirling's formula one finds, in the notation above, 
i/j(x) = x a , ^{x) ~ <j>*{x) ~ x 1/a , 

and so 

- logP(X > x) ~ <j>{x) ~ (1 - a)a Q/(1 " a) a; 1/(1 - a) . 

This result is derived twice in [6], p. 350, 352, and again in [5], Th. 8.1.12, 
as an instance of Kasahara's Tauberian theorem. For the mgf itself, one has 

ML a (t) ~exp{t 1/Q }/a 

(Cartwright [8], p. 50). The Mittag-Leffler law of index a is closely connected 
to the stable subordinator (increasing stable process) of index a, and the 
result above translates into a tail estimate at zero for X — X(l) in the stable 
subordinator: 

-logP(X < x) ~ (1 - a)a a/ ^~ a) /x a/ ^- a) (x | 0); 

see [6] 349-350, [5] Section 8.3.4. 

The limiting case p — 1 relates to entire cfs of infinite order. This is 
a boundary case of the above, and may be handled by using second-order 
theory of regular variation (de Haan theory: [5] Chapter 3). Let be a 
positive function with <p(t)/t — > oo as t — > oo, and similarly for a. Suppose 
the second-order condition 

6(tx) — x6(t) , . 
VK / \ ^xlogx (t^oo) Vx>0 
ta(t) 6 v ; 
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holds. Then if the tail function is everywhere positive and the cf is entire, 
the following are equivalent: 



lim inf 



logP(|X| > x)-<t>{x) 
xa(x) 



= 



for some 0, a as above, 



lim inf 



[E(\xn\ 




= 



n^oo 



a{n) 



for some ip and some a G II (the de Haan class II; see e.g. [5] Ch. 3). Then 



and a ~ a. We refer to Geluk [15], Geluk, de Haan and Stadtmiiller [16] for 
details. 

3. N on- analytic characteristic function 

A cf is defined in the first instance on a line in the complex plane (the real 
axis if one includes the T in the definition, the imaginary axis if, as here, one 
does not), but can be continued analytically to a strip in the complex plane 
(which may be the whole plane, a half-plane, a strip of positive finite width 
or the imaginary axis). A cf is analytic iff all moments fi n are finite, and 
limsup n _> 00 [|/i n | / M!] 1 / n = 1/R is finite (R is then the radius of convergence 
of the Taylor series of the cf, or mgf). Equivalently, it is analytic iff for some 
R > the tail function T(x) := P{\X\ > x) is 0(e~ rx ) for all r < R, and 
then the strip of regularity contains the strip | Im z\ < R ([29], Th. 7.2.1). 
We proceed to consider results analogous to those above, but with the tails 
heavy enough to cause non-analyticity - that is, with the log-tail decaying 
more slowly than a multiple of x. 

A recent result of this kind, which may serve as a prototype, is due to 
Konig and Morters ([23], Lemma 2.3): for X non-negative, p > 1 and c real, 



with a similar implication with limsup in place of limit. (In their setting, 
p = 1,2, . . ., and they consider the intersection of p independent Brownian 




n^oo n 




implies 



lim x~ 1/p \ogP(X > x) 
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motions in d dimensions. We exclude here the case p — 1, which falls under 
analyticity above, and is not relevant to their intersection problem, but there 
is no need here to restrict to p integer.) Their proof is simple, using only 
Markov's inequality, Stirling's formula and reduction of the limit to a suitable 
sequence. On the other hand, the setting of their main results makes heavy 
use of large-deviation theory. 
4- Complements 

1. Laplace's method. The essence of Tauberian theorems of exponential 
type, and of large deviations, is that growth is dominated by behaviour in 
the neighbourhood of the maximum. This is simply a reflection of how over- 
whelmingly rapid exponential growth (or decay) is. The prototype in analysis 
is Laplace's method (for which see e.g. de Bruijn [7]), which goes back to 
Laplace's Theorie analytique des probabilites of 1820 and of which the most 
telling instance is Stirling's formula of 1730. Its applicability in statistical 
mechanics rests ultimately on the enormous magnitude of Avogadro's con- 
stant, c. 6 x 10 23 . 

2. The Wiman- Valiron method. In the discrete analogue of Laplace's method 
(as in [6], 346), one examines the growth of an entire function (an entire cf in 
our case) by localizing at the term of maximum modulus in its Taylor series. 
For surveys of the Wiman- Valiron method, see Hayman [18], Fenton [13]. 

3. Esscher transforms. For F a law with mgf M(t), the law F t {dx) : = 
e tx dF(x)/M(t) is called the associated law by Feller [12], XVI. 7, or the ex- 
ponentially tilted law. Its use goes back to F. Esscher in 1932, and later 
to Cramer in 1938; it is accordingly called the Esscher transform. It is 
widely used in ruin theory and insurance and actuarial mathematics; see e.g. 
Mikosch [31], Section 4.2. In brief, one may use it to change at will the point 
in whose neighbourhood one wishes to probe the asymptotic behaviour. 

4. Mellin transforms. As mentioned earlier, with all moments finite as here, 
one may consider the Mellin transform J-(s) := E[X S ] of X, and seek to link 
its behaviour for large s with the tail behaviour of X. One early example 
is the work of Wagner [39], also considered in [6] (the beginning of my col- 
laboration with Jef Teugels, which led to [5], the origin of Fenchel duality 
in regular variation, and one source, together with [1], of smooth variation - 
see [5] Section 1.8). 

5. Branching processes. Tauberian theorems of exponential type are also 
useful in branching processes; see e.g. [2]. Again, the theme is that when 
populations grow, they grow so rapidly that the contribution near the max- 
imum is all that counts. 
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6. Random potentials. We close by mentioning the wealth of related liter- 
ature on random Schodinger operators in models of disordered media, etc. 
For a recent survey, with many references, see Leschke et al. [27]. It is no 
accident that Tauberian theorems of exponential type were also discovered 
in the theoretical physics literature, as the Minlos-Povsner Tauberian theo- 
rem [32], [33]. As Cindy says, all good mathematicians are also interested in 
physics. 

Postscript. It is a pleasure to record a personal debt to Cindy Greenwood 
here. The long Chapter 4 in [5] (64 pages), entitled Abelian and Tauberian 
theorems, originally covered Mercerian theorems also. Only after repeated 
and urgent pleas from Cindy did I detach the Mercerian component, which 
became Chapter 5 (25 pages - with Chapter V in Pitt [35] , incidentally, the 
only textbook chapters I know, though see also Korevaar [24] IV. 13, 14, Paley 
and Wiener [34] IV. 18). She pointed out with remorseless clarity that sepa- 
rated, the Mercerian material would be impossible to miss, but unseparated 
risked sinking without trace. I wish I had had the courage of her convictions 
and split the Abelian and Tauberian parts also. 

Acknowledgement. I thank Peter Morters for rekindling my interest in this 
subject. 
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